
figure
hold on
%x = linspace(-1, 1, 100); % 生成0到2*pi之间的100个点
plot(uhG(:,3),'r*');
plot(ureal(:,3),'k-');
legend('MRWeno限制器','解析解');
hold off
